Quantum theory of frequency pulling in the cavity-QED microlaser 
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The spectrum of the cavity-QED microlaser/micromaser is expected to show distinctive features 
of the coherent light-matter interaction, which are obscured in the conventional Schawlow-Townes 
linewidth theory. However, the spectral studies has been limited to resonant atom-cavity interaction 
so far. Here we consider the dispersive interaction in the off-resonance case, from which we uncover 
a quantum frequency pulling effect in the microlaser/micromaser spectrum. We present a quantum 
theory of the spectrum which introduces the notion of a frequency-pulling distribution associated 
with the photon number. In contrast to the conventional laser, periodic variation of the mean 
frequency pulling is observed with increasing pump parameter and it is attributed to the strong 
atom-cavity coupling. The pulling distribution gives rise to a spectral broadening, which can be 
dominant over the non-dispersive broadening addressed in the previous works. We also developed 
a corresponding semiclassical theory and discuss how the introduction of the frequency shift fits in 
with the extended quantum theory. 

PACS numbers: 42.50Ar, 42.50Nn, 42.50Pq, 42.55Ah 
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I. INTRODUCTION 

The microlaser/micromaser [HQ operates with a small 
number of atoms under the strong atom-cavity cou- 
pling condition of the cavity quantum electrodynamics 
(QED). The manifestation of atom-field coherent inter- 
action, which is masked by stochastic averages in the 
conventional laser, have been observed in the micro- 
maser/microlaser [3|— l5|j . Moreover, its quantum particle 
property has revealed the nonclassical aspect of the elec- 
tromagnetic field . The wave property, on the other 
hand, which is expected to exhibit the complementary 
nature still awaits its experimental verification. In par- 
ticular, the power spectral density or the first-order co- 
herence is expected to show features [9] unexplainable 
by the conventional Schawlow-Townes linewidth theory 
[lOj . Since the seminal work on the theory of the mi- 
cromaser spectrum Q, more elaborate and refined cal- 
culations [llHlq ] have followed. Indirect measurement of 
the micromaser spectrum has been proposed by using an 
atom probe |16l4l8l |. In the microlaser operating in the 
optical region, the spectrum measurement is expected to 
be rather straightforward. 

It should be noted that all those previous works on the 
micromaser/microlaser spectrum have been restricted to 
the atom-cavity resonance condition, and thus they can- 
not be applied to off-resonance cases by simple modifi- 
cations because dispersive interactions can induce a fre- 
quency shift, resulting in a qualitatively different spectral 
lineshape. Even in the conventional laser, for example, 
the phenomena called frequency pulling or mode pulling 
makes the true laser oscillation frequency shifted from 
a passive cavity resonance frequency in the presence of 
a frequency detuning between the gain medium and the 
cavity. It is thus natural to ask how the spectrum changes 
as we remove the constraint of resonance over the whole 
tuning range of the microlaser. 

Here we present an extended quantum theory of the 



microlaser spectrum including off-resonance cases. Our 
theory predicts a quantum frequency pulling effect, a 
novel feature not found in the conventional laser. More- 
over, the amount of the frequency pulling strongly de- 
pends on the pumping parameter owing to the strong 
atom-cavity coupling. In the standard laser theory of 
Lamb and Scully (19J, in contrast, the amount of the 
pulling is solely determined by the linewidths of the 
atomic and the cavity resonances independent of a pump- 
ing parameter. In particular, the unitary evolution of the 
atomic state gives rise to a periodic characteristic in the 
pulling. The analysis given here is also applicable to 
the microwave counterpart as long as thermal microwave 
photons are negligible. 

The quantum frequency pulling to be discussed be- 
low is indexed by the photon number, the uncertainty of 
which leads to spreading in the frequency shift or a spec- 
tral broadening. This mechanism of dispersive broad- 
ening has never been addressed so far although it can 
easily exceed the nondispersive contribution in the spec- 
tral width under the present experimental conditions of 
the microlaser. Only microscopic lasers like the micro- 
laser/micromaser are adequate to study this quantum 
dispersive effect: in the conventional laser, the differen- 
tial pulling by one- photon increment/decrement is neg- 
ligible and thus not a noticeable effect. We note that 
laser oscillation frequency shift with respect to the cav- 
ity has been discussed in the socalled mazer [2(| , a maser 
pumped b y th e matter wave, and in the trapping state 
spectrum [111 ], but not in a general context as in the 
present paper. To help gain an insight on the frequency 
pulling effect, we have also present a corresponding semi- 
classical theory. 

This paper is composed of the quantum theory in Sec. 
HI! and the semiclassical theory in Sec. IIII1 The quantum 
theory is based on the quantum regression (Sec. Ill A"| and 
the master equation (Sec. Ill Bj) . Central analytic formulae 
are derived in Sec. Ill Cl and Appendix |A"1 After examining 
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the physical meaning of the result of the quantum theory 
and its conventional laser limit (Sec. IIIDI) . typical spec- 
tral lineshapes (Sec lII E| ) and the dependence on detuning 
fSec lIIF[) and pumping fSec lII Gl) arc investigated. The 
semiclassical theory presented in Sec. IIII Al is based on 
the Maxwell-Schrodinger equation and is compared with 
the quantum theory in Sec. IIII Bl The correction made 
by a newly-introduced frequency variable in the semiclas- 
sical theory is elucidated by using a graphical method in 
Appendix [B] 

II. QUANTUM THEORY 

In this section we generalize the master-equation-based 
theory of the micromaser spectrum Q in order to include 
the effect of atom-cavity detuning. The key expressions 
relating steady-state variables are derived analytically. 
The numerical solution of the master equation is used 
to determine these steady-state variables and, in turn, 
to perform exact calculation of the quantities related to 
the frequency pulling. The quantum distribution of the 
frequency pulling is naturally introduced along this line 
of analysis. 



Therefore, (a'(t)a(0)) is the expectation value of the 
field amplitude subject to the transformed density ma- 
trix p(t) = tt R [U{t)a{Q)x(S>)U\t)]. According to the 
quantum regression theorem [19| , the equation of motion 
for p(t) is given by that of the pit), that is, 

m = cm- (6) 

In the photon number basis the first-order coherence 
function is related to the off-diagonal element of the den- 
sity matrix p n ,n+i = P~n as 

( a t(t) (o))=x;v^+T^ i) (*) s (7) 

n 

where pn\t) is the solution of Eq. ([5]). The initial value 
of the transformed density matrix is given by 

p(0) = tr R [a(0) X (0)] = a(0)p(0), (8) 

which is in the photon number basis 

(n\p(0)\n + 1) = Vn+l(n + l|p(0)|n + 1) (9) 

with all other terms being zero. 



First-order correlation function 



B. Master equation for p\ 



(i) 



The spectrum of an optical field is obtained by taking 
the Fourier transform of its first-order correlation func- 
tion which is written as 



(a)(t)a{Q)) = tr FeJt [ X (0)ot(t)a(0)] , 



(1) 



where x is the density operator of the field(F)- 
reservoir(R) system and and a are the creation and an- 
nihilation operators, respectively, of the laser field. The 
time origin t = here is set to an arbitrary moment at 
which the steady state of the system is achieved. We 
have knowledge on the evolution of the field subsystem, 
which is described by 



P {t) = tlR [ X (t)} 

via the master equation 

p(t)=Cp{t), 



(2) 



(3) 



the detailed form of which will be described below. The 
time evolution of the total system is given by U (t) such 
that at(i) = [7t(t) a t(0)f7(t). We proceed as 

(ot(t) ffl (0)) = ti mR [ X (0)U\t)a^0)U(t)a(0)] 
= tr F [a+(0)tr fl [U(t)a(0) x (0)U^t)]} 
= tr F [ot(0)p(t)], (4) 



to give 



The next task is to find the time evolution of pn ^ — 
Pn,n+i- The master equation is written by considering 
the change of the field due to the coherent interaction 
with the initially inverted atoms and the incoherent de- 
cay. We assume the tophat interaction model in which 
the atom-cavity interaction parameterized by the con- 
stant atom-cavity coupling g persists for a fixed amount 
of time t for all atoms that traverse the cavity mode. 
The mean number of atoms inside the interaction region 
is designated by N. The realistic modification such as the 
spatial variation of the cavity mode and the inhomoge- 
neous distribution of the interaction time can be readily 
included in the numerical calculation. 

The coarse-grained master equation for the microlaser, 
written in the photon number basis, reads [l9| 

Pnra Q>nmPn,m ^nmPn+l,m+l ~l~ CnmPn — l,m — 1 (10) 

where 



^nm — ' a COS [ — z — I COS 



A' 2 . fn n T\ . fn mT 

sm — — sin 



2 



A' fn n r 

, .A' /fi m r\ . fn n r 



(at(t)a(O)) = (J) m . 



(5) 



7 c (n + to), 

(11) 
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b nm = 2 7c V(n + l)(m+l), (12) 

(13) 

with the /leM-atom detuning A' = to — luq, the n-photon 
Rabi frequency O n = ^4g 2 (n + 1) + A' 2 , the injection 
rate of the atomic beam r a — N/t, and the cavity de- 
cay rate j c (a half width). Note that the frequency of 
the laser field uj does not necessarily coincide with the 
passive cavity frequency lu c if we are open to possible 
frequency shift. In that sense we reserve A for the pas- 
sive cavity detuning A = lj c — ujq . Since we are interested 
in the steady-state operation, in which the frequency of 
the field is already shifted, the master equation is writ- 



ten in terms of the resulting shifted detuning A'. The 
introduction of A' and the corresponding modification of 
fi„ in Eqs. (fTT]) - (fT5]) are the key changes made to the 
quantum theory of Ref. 9]. The imaginary part of Eq. 
(|10p. which appears only in the presence of a nonzero 
detuning, already suggests the introduction of frequency 
pulling since the steady state p nm = cannot be achieved 
without the frequency shifting transformation such as 

P'nm(t) = Pnm(t)e^- m ^-^. 

C. Complex eigenvalue of pn and the emergence 
of frequency pulling 

The equation of motion for the off-diagonal element 

(i) 

Pn appears as 



-• W £^'," + " *- (^) *»> (^) & + ^ v^wtt^ u (») 



where we defined 



fl n r\ ( Sl n+ ir\ A' 2 . / fi„r\ . ( O n+ ir\ . A' / £1„t\ . /Q n+1 T 



2^ 1 I 2 J I 2 J fi n f2 n _i_i V 2 I \ 2 I f2«+i 



cos ( Qn+lT \ sin f°" T> ) + 4g 2 V(n + l)(n + 2) ^ fQ n T 

fin \ 2 / \ 2 J 



(^)sin(^)-l} + 2 7c y^TT) 
7c(2n + l), (15) 



and we rearranged the terms in order to employ the 
detailed balance between 2"f c y/ (n + l)(n + 2)p| l 1 | 1 and 

r ° 4g X ( „"n„ 1 | ( 1 n+2) sin (V) sin (^T 1 ) p£ } in the stead y 

state Such a balance guides pn to decay exponen- 
tially with a single time constant \p n . The detailed bal- 
ance holds for most of the parameter space except the 
trapping state condition 0, [2l[ , at which the truncation 
of the Hilbert state occurs. We note that the single-decay 
ansatz is known to give more or less an inaccurate value 
of the linewidth for the trapping state [HI , and the cor- 
rect spectrum often comes up with a nontrivial lineshape 
[TT| . We exclude this singular case in our discussion for 
simplicity. As noted in the definition of A' and A in Sec. 
Ill Bt A' = A + (6 n ) is understood in the steady state. 



The corresponding field correlation function is then ob- 
tained from Eq. ([7]) with t > as 

(at(f)a(Q)) = ^v^TT^We-^" 4 

n 

= ^(n + lWie-W 
n 

Egn + lWe-^' (16) 

n 

where D n and — 5 n are defined as the real and imagi- 
nary parts of \p n , respectively, and p n is the steady- 
state photon number distribution which can be obtained 
by the standard theory of the micromaser [g^]. The op- 
tical spectrum obtained by 

1 r°° 

S(v-to c ) = -Re (aUt)a(0))e- i ^- u ^ t dt (17) 

I" JO 

is therefore the sum of individual Lorentzians centered at 
v = co c + 6 n with a half linewidth D n . Here D n is iden- 
tified with a phase diffusion constant at the fixed oscilla- 
tion frequency while 6 n is interpreted as a frequency shift. 
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We note both D n and 5 n depend on n, and thereby have 
distributions associated with the photon number. While 
the phase diffusion D n represents a spectral broadening 
by itself, the frequency shift S n gives rise to an additional 
broadening by spreading out the center frequencies of the 
individual Lorentzians. Note that the contribution from 
each Lorentzian is weighted by np n . 

In the limit of n ^> 1, the real part D n is further 
simplified as 



D n ~ 2r sin 



7c 
An 



(18) 



as shown in Appendix [X] This result appears as a simple 
generalization of the on-resonance result [j| by replacing 
2g\/n + 1 in the denominator of the sine term with tt n . 

Let us examine the imaginary part S n which has never 
been addressed so far in the literature. In the same limit 
of n ^ 1, it is simplified as (see Appendix IA")l 



1 ni 



i - 



sin(J7 n r) 

fl„T 



(19) 



Notice that, on resonance (A = 0), the frequency pulling 
8 n should vanish to satisfy Eq. (fT9"l) . The frequency 
pulling is truly a result of dispersive atom-cavity inter- 
action, which induces an additional optical path length 
other than that of the cold cavity. The overall minus sign 
indicates that the sign of S n is opposite to that of the de- 
tuning so that the oscillation frequency of the microlaser 
is pulled towards the atomic resonance. 

Since fi„ in the right-hand-side of Eq. (fT9"|) also con- 
tains (S n ) in the effective detuning, Eq. (fT9]) can not be 
written in a closed form equation for S n . Although some 
of its properties can be inferred from an analytic form as 
we will do below, the exact calculation of the frequency 
pulling is presented by numerical means mostly with the 
arameters relevant to the microlaser experiment of Refs. 
"J throughout this paper. In the numerical calculation 
the equivalent Gaussian mode function is employed to be 
more realistic. 



D. Physical meaning and the conventional laser 
limit of the frequency pulling 

Before we proceed with systematic numerical calcula- 
tions, let us appreciate the physical meaning of Eq. (|T9|) . 
The mean value of the frequency pulling is approximately 
equal to d n with n ~ (n) . We can rearrange the terms in 
Eq. OSJ) as 



iy A' N 

<^} = -^rx-x^) 



(20) 



where the characteristic function £($) (plotted in Fig.QJ 
is defined as 
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FIG. 1: The characteristic function representing the 

Rabi angle dependence of the frequency pulling. 



Here we 
27c/(.9 2 t 



the threshold atom number Nth 
the semiclassical Rabi angle $ 
n) + A' 2 t = fir assuming (n) ^> 1. 



used 
and 



(21) 



The first multiplicative factor in the right-hand-side of 
Eq. ([20|) corresponds to ■y r A/2^„, the frequency pulling 
in the conventional laser [19j , if we identify the dephasing 
rate of the medium j p with the transit-time bandwidth 
(2r) _1 . It has a constant value determined by the rel- 
evant bandwidths of the cavity and the gain medium, 
independent of the pumping rate. In the case of the mi- 
crolaser we have additional multiplicative factors in Eq. 
(HO]). 

The factor N/N t h, equivalent to the effective pumping 
parameter often denoted by 6 2 = Ng 2 r/2"/ c |22], shows 
the increase of the frequency shift by pumping while £($) 
term indicates the reduction resulting from a Rabi an- 
gle associated with the growing laser intensity. In the 
conventional laser these opposite tendencies balance each 
other to result in a constant amount of the pulling inde- 
pendent of the pumping. In other words, the gain sat- 
uration approximately leads to N oc (n) cx <I> 2 / 'g 2 above 
the threshold while £($) ~ makes N£($) constant. 

In the microlaser, on the other hand, the gain is os- 
cillatory with (n) so that we encounter the reduction of 
gain well above the threshold and the corresponding sat- 
uration of (n) as a function of ./V [see Fig. fTUT a) below]. 
The pumping by the atomic flux N is thus translated into 
the increase of the amount of pulling. 

The other pumping channel, the coupling constant g, 
brings the mean photon number down when it is in- 
creased, and the corresponding variation of $ 2 with re- 
spect to g is rather small. Thus g 2 behavior contained in 
N/Nth term is dominant also in this case. The quadratic 
dependence comes from the feedback mechanism inher- 
ent in the laser amplification process, in which the atomic 
polarization induced by the cavity field via coupling g 
drives the cavity field again by the same coupling. 
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The notable exception of the general tendency of Ng 2 - 
type increase in the frequency pulling happens at the 
socalled multiple thresholds, which are associated with 
the periodic Rabi flopping [4{. In this case, (n) changes 
discontinuously, and so do $ and resulting in the 

periodic rise and fall of the pulling with 8 2 . Thus the 
microlaser operation is segmented into branch solutions 
associated with multiple thresholds. 

The discontinuous change of the system can be ex- 
plained in terms of the coherent evolution of the atomic 
state. As the system approaches the saturation regime 
out of the linear regime of each branch solution, the pre- 
cession of the Bloch vector representing the internal state 
of the gain medium reduces the emission probability of 
photon from each atom while the phase shift per atom 
grows. At the transition to a higher branch solution, the 
system is brought to a new steady state by making the 
precession angle of the Bloch vector exceed 2ir. In the 
middle of such interaction the atomic state traces back to 
its initial state and starts the evolution all over again. In 
doing so, the emission probability is boosted up and re- 
sults in a sudden rise of the mean photon number - also 
known as a quantum jump |3[ - while the phase shift 
per atom diminishes. Hence, upon the multiple thresh- 
olds, the microlaser behaves in a counter-intuitive man- 
ner in the sense that it becomes less refractive with a 
more dense medium of atoms. This nontrivial pumping 
dependence based on the strong atom-cavity coupling is 
one of the distinctive features of the microlaser. The 
argument given here will be examined in detail by nu- 
merical calculations in Sec. Ill Gl 

The result of the conventional laser can be recov- 
ered if we average the microlaser result over a ficti- 
tious dwell time distribution given by P(t) = 2 7p e _27pT 
which spoils the well-defined atom-cavity cou- 
pling by introducing stochastic interaction time. Us- 
ing J™ re- 2 ^ (l - ^) dr = and assum- 

ing the pulling to be small compared to the detuning 
(A' ~ A) we have 



(Sr. 



g 2 A 



'2 7p (4 7p 2 + 2 ) 2 7p 



(22) 



where the last equality employs the steady state condi- 
tion of the conventional laser 4 ^ r 2 + n 2 = 2 7c . 



E. Spectral lineshape 

The overall lineshape is given by a sum of Lorentzian 
curves spread by the photon number distribution np n . 
The amount of the spread 6 n as well as p n depends on 
the detuning A, thus we expect a significant change of 
the lineshape as we tune the cavity. The overall spectral 
shift from the passive cavity resonance is given by S n with 
n corresponding to a predominant p n . To obtain the de- 
tailed spectral distribution we have to take into account 
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FIG. 2: (a) The photon number distribution p n (solid curve), 
the corresponding D n (dotted curve) and S n (dashed curve) 
calculated with gr = 0.124 and 7c r — 0.049, corresponding 
to the experimental parameters of Refs. 0, ||- The tuning 
parameters are fixed at N = 100 and A/7 C = 14.76. (b) The 
normalized spectrum as a result of superposing the Lorentzian 
curves corresponding to different n's. The major spectral en- 
ergy is confined in a frequency range lower than v = cj c . The 
constituent Lorentzian curves with An = ±10 are depicted in 
dashed curves for comparison. The spectrum is normalized 
to a unity height. 



both the intrinsic linewidth D n and the distribution of 
the shift S n - 

Let us for now assume a single-peaked p n as is the case 
for most of the microlaser parameters. To a good ap- 
proximation, p n with a large (n) is well described by the 
Gaussian curve. The constituent Lorentzian curves have 
little variation in their widths D n within the range of 
substantial probability p n as in Fig. [2ja) . Thus the spec- 
tral lineshape is approximately given by a Voigt integral 
centered around lj c + (6 n ) as in Fig. Hfb). Of course the 
lineshape exactly reduces to a Lorentzian on resonance. 
The pulling-induced broadening exceeds the linewidth of 
any constituent Lorentzian even with a detuning of sev- 
eral times of 7c away from the atomic resonance, thus 
becomes the dominant source of decoherence for a large 
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part of the tuning range of the microlaser. 

One may be tempted to use this n-dependent shift 
to obtain information on the photon number distribu- 
tion, which is possible if the differential pulling is much 
larger than the individual linewidth D n . Unfortunately, 
we could not identify a parameter set satisfying this con- 
dition because any parameter sets we tried not only in- 
creased the differential pulling but also the individual 
linewidth as well. 

However, the photon number distribution can do a sig- 
nificant role in the spectrum at some particular operat- 
ing points. It is well known that the photon number 
distribution is not single peaked if the system under- 
goes any of the multiple thresholds as shown in Fig. [3ja) 
[3,0]. Although S n is a slowly varying function of n, the 
large separation between the two possible solutions of 
p n is enough to separate the spectral peaks as depicted 
in Fig[3]Jb). In this case, the narrower photon number 
distribution in the lower branch solution near n ~ 500 
corresponds to a broader spectrum near v — uj c ~ —4.25. 
This is because the higher branch solution, despite the 
relatively broad p n , generally corresponds to an order- 
of-magnitude sharper spectral peak D n as well as much 
smaller 8 n . In other words, the spectral energy is redis- 
tributed according to S n and D n while the ratio of the 
integrated area under each peak of np n , about 1.2:1 as 
shown in the inset of Fig. [U[b) , is preserved in the spec- 
tral domain. In reverse, one should also be able to obtain 
the binary p n distribution from the observed dual-peak 
spectrum as long as D n and 5 n in each peak of p n are 
approximately constant. 
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F. Detuning dependence 

The detuning dependence of the spectrum confirms 
that the oscillation frequency of the microlaser is indeed 
pulled toward the atomic line center as in the conven- 
tional laser. The magnitude of the effective detuning A' 
is always smaller than that of the passive detuning A. 
The spectra calculated with three different detunings are 
depicted in Fig. [4] with respect to the passive cavity fre- 
quency. We note that the center frequency of the spectra 
moves the farther away from lo c for the larger detuning. 
The lineshape transforms from a Lorentzian-like curve to 
a Gaussian-like curve as the frequency excursion induced 
by 5 n gets larger than the individual width D n . 

Let us now examine the mean (S n ) and the variance 
AS n of the pulling distribution systematically. Before 
the microlaser reaches the second threshold, i.e., within 
the first branch solution, the mean spectral shift (S n ) 
as a function of A follows a familiar dispersion curve 
[Fig. Ufa)]. The spectral broadening AS n is accordingly 
maximized at the detuning which maximizes (5 n ) [Fig. 
Hfb)]. There we used the standard deviation of S n as a 
measure of the linewidth, that is, A8 n = \J (8%) — (S n ) 2 . 
We note that the relative contribution of D n becomes 
significant on resonance, for which the contribution of 8 n 



FIG. 3: Plot of the same quantities as in Fig. [2] except for 
the tuning parameters set to a multiple-threshold transition 
point, N = 367 and A/*y c = 16.61. The inset in (b) shows 
the weight function (n + l)p n +i. 
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FIG. 4: The spectral lineshape depending on the detuning 
A/7 C =11.87 (solid), 23.88 (dashed), 48.29 (dotted). The 
system parameters are gr — 0.124, 7 c r = 0.049, and N = 100. 
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FIG. 5: (a) The mean photon number (n) (dashed) and the 
mean frequency pulling {8„)/"/c (solid), (b) The spectral 
broadening (D n )/^y c (dotted) and A8 n /"/ c (solid) as the pas- 
sive cavity frequency is tuned to near resonance. Plotted with 
gr = 0.124, -y c r = 0.049, and N = 10. 



FIG. 6: Same as Fig. [5] except for N = 300 to reach the 
saturation regime. 



tive explanation previously given in Sec. Ill Dl 



just vanishes. The contribution by D n is also dominant at 
far-off resonance, which corresponds to a below-threshold 
regime; since the field is just the atomic emission, not 
a lasing field with feedback, filtered by the cavity, we 
confirm (D n ) ~ j c there. 

In the saturation regime within each branch solution, 
on the other hand, the curve (S n ) vs. A turns into a 
power-broadened dispersion curve (Fig. [5]) . The contri- 
bution of the dispersive broadening becomes dominant 
almost over the whole tuning range of the microlaser. 
Thus the spectrum cannot be explained without the no- 
tion of the frequency pulling. 

An extraordinary dispersion curve can be obtained if 
we pump the system (see also Sec. Ill Ol below) harder 
so that more than one branch solution can be attained 
within the cavity-tuning range. The mean photon num- 
ber or the output intensity of the microlaser as a function 
of the detuning is known to form a set of discontinuous 
fragments of individual branch solutions 0, HH as shown 
in Fig. [7{a). At the very point of the transition where 
the fragments are connected, the system is bistable as 
discussed above in Fig. [3] Note that the higher branch, 
which occurs at the smaller detuning, has the larger mean 
photon number. The detuning curve of (S n ) is also a dis- 
continuous stepwise dispersion curve. The amount of the 
pulling (8 n ) is, however, smaller for the higher branch in 
Fig. 13a). This opposed behavior of (n) and (S n ) upon 
the multiple thresholds is in agreement with the qualita- 



G. Pumping dependence 

As we increase the pumping at a fixed detuning A, 
the effective detuning A' tend to diminish monotonically. 
However, it is interrupted by the aforementioned reset- 
ting of the atomic state, which is subject to the periodic 
Rabi flopping. Like other characteristics of the micro- 
laser such as (n) or Mandel Q, the frequency pulling also 
undergoes a periodic change as a function of pumping. 

The pumping parameter of the micromaser is given by 
= \/N ex gT where N ex — N/2j c t. As discussed in Sec. 
Ill Dt we expect the overall pumping dependence to be 
9 2 and its periodic recurrence at the multiple thresholds. 
The pumping dependence in the standard micromaser 
theory [2l|, [22[ is usually presented as gr is scanned at 
a fixed atomic flux whereas scanning N ex at fixed gr is 
more relevant to the experimental practice of the micro- 
laser gg. 

1. Increasing the atom-cavity coupling 

We tune the coupling constant g at a fixed atom num- 
ber N in Fig. El The multiple thresholds are found for all 
depicted variables. The mean photon number maintains 
its typical feature [22j regardless of the correction by the 
presence of 8 n . The increase of the frequency pulling 
in the first branch solution is mildly nonlinear in Fig. 
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FIG. 7: Same as Fig.[5]except for N = 700 to reach the second 
branch solution near resonance. 
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FIG. 9: The pumping dependence in terms of the atom-cavity 
coupling constant g at the fixed detuning A/"/ c = —15.0 and 
N = 250. (a) The mean photon number, (b) the mean fre- 
quency pulling and (c) the spectral broadening {D„)/7 C (dot- 
ted) and A(5„/7 C (solid). The arrow indicates the magnitude 
of the maximum achievable g in the current experimental set- 
ting of the microlaser [8j. 



it drops to a smaller value. Up to the first branch solu- 
tion the dispersive broadening dominates over the non- 
dispersive one. 



FIG. 8: Same as Fig. [5]except for N = 1300 to reach the third 
branch solution near resonance. 

Mb) due to the aforementioned g 2 effect. The curvature 
slightly changes at the first threshold of gr ~ 0.02. The 
amount of the pulling induced by a small number of gain 
atoms is enhanced by the large coupling constant. The 
spectral broadening generally becomes larger for stronger 
pumping except in the multiple-threshold regions, where 



2. Increasing the atomic flux 

The pumping by the atomic flux also brings about the 
multiple thresholds as shown in Fig. [TO] The mean pho- 
ton number is stepwise as verified in the previous experi- 
ments dLH, [HI- The frequency pulling also increases here 
within each branch solution. It keeps increasing overall 
except for the periodic drops. The spectral broadening 
is dominated by the frequency pulling for the parameters 
shown in Fig. [TU] 
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A. Maxwell-Schrodinger equation 
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FIG. 10: The dependence on the pumping via the atomic flux 
(mean atom number N) at the fixed detuning A/7 C = —15.0 
and gr = 0.124. (a) The mean photon number, (b) the mean 
frequency pulling and (c) the spectral broadening (D n )/~f c 
(dotted) and A<5 n /7 C (solid). 



III. SEMICLASSICAL THEORY 



A set of coupled equations describing the microlaser 
is comprised of the Maxwell equation for the cavity 
field driven by the induced atomic polarization and the 
Schrodinger equation for the unitary evolution of the 
atomic state under the electric dipole interaction. In 
the slowly varying envelope approximation the Maxwell 
equation reduces to 



£(t)+j c £(t) =i 



2t7LU c 
V 



P(t' 



dt' 



(23) 



and the Schrodinger equation to 



T{t) - i(u c - wo)P(t) = -»^£(t)r(t), (24) 



and 



lm[£*{t)V(t)} 



(25) 



where £ and V are the slowly varying amplitudes of the 
cavity electric field and the induced atomic polarization, 
respectively, V is the mode volume of the cavity field, fi 
is the induced atomi c dip ole moment and r is the atomic 
population inversion |24l25| . Whereas the previous stud- 
ies left the variables £ and V complex [24], [25| , we find it 
convenient to decompose them into an amplitude-phase 
form to highlight the role of the in-phase and quadrature 
components of the induced dipole. With normalization 
in mind we rewrite 



£{t) 



V 



■a(t)e' 



i<t>{t) 



V(t) = t i[S 1 (t)-tS 2 (t)]e-^ 



(26) 



(27) 



where a(t), 4>{t), S\{t) and S%(t) are all real. We rename 
r(t) = Ss(t) in the context of the Bloch vector. The 
corresponding Maxwell equation then reads 

d(t) + 7o o(t) = Ng j T ^P-dif, (28) 



In this section, we construct a semiclassical theory for 
intuitive understanding of the frequency pulling phenom- 
ena. The source term of the Maxwell equation clarifies 
the origin of the frequency shift in terms of the Bloch vec- 
tor. While the semiclassical theory of the conventional 
laser is comprised of a pair of equations for both inten- 
sity and oscillation frequency, the semiclassical equations 
concerning the microlaser have so far been limited to the 
intensity alone @, [H[ . Here we provide the other equa- 
tion for the frequency shift as well as consequential modi- 
fication of the intensity equation. The correspondence to 
the quantum theory is also made in the mean frequency 
pulling. 



* {t) = - Ns [Wrr d,< - (29) 

and the Schrodinger equation turns into 

S 1 (t)-(A + j>)s 2 (t)=0, (30) 

S 2 {t) + (A + ^S 1 (t)=2gaS 3 (t), (31) 

S 3 (t) = -2gaS 2 (t), (32) 
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where the atom-cavity coupling constant is given by g = 
- :r "" ,; Here the time scale involved is different for the 



v 



field variable and the atomic variables; the atomic state is 
transient during a single transit through the cavity mode 
while the field evolves in a longer time scale over many 
atomic transit events. We observe that the amplitude of 
the cavity field is driven by the quadrature component S2 
of the atomic dipole while the frequency shift is produced 
from the in-phase atomic dipole Si . The equations for the 
atomic variables are exactly the optical Bloch equation 
except for the shifted field frequency A' = A + (j>. 

With the values of a and 5 = 4> held constant in the 
steady state, we can readily find the solution of the equa- 
tions of the atomic variables, Eqs. (|30" l) -(f32 j) . The result 
of this solution will be used to eliminate the atomic vari- 
ables in Eqs. (|28f) and (l29t below, arriving at the coupled 
equation of a and 5 for the determination of their steady 
state values. 

With the initial conditions which represent the perfect 
population inversion of the incident atoms, 53(0) = 1 
and 5*1(0) = 52 (0) = 0, the set of solutions is found as 



Slit) 



2ga(A + 5) 



Ag 2 a 2 
X [l - 



(A + 6) 2 

os (VVa 2 + (A + <5) 2 t)] , (33) 
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FIG. 11: The comparison of the quantum (thick curve) and 
the semiclassical (thin curve) calculations for (a) the mean 
photon number and (b) the mean frequency pulling. The 
same parameters are used as in Fig. [7] 



s 2 (t) 



and 



Ss(t) = 



2ga 



y/4g 2 a 2 + (A + 5) 2 



sin ( VVa 2 + (A + S) 2 tj , 
(34) 



1 



4 5 2 a 2 + (A + 5) 2 

x [Ag 2 a 2 cos (VVa 2 + (A 4- 5) 2 t 



(A + 5) 2 }. 
(35) 



B. Equations for intensity and frequency 



The integration of Eq. (|28|) in the absence of S pro- 
vides the conventionally known gain- loss equation [4ll23l]. 
In the generalized approach given here, the dispersive 
correction to the microlaser gain is made via the effec- 
tive detuning A'. More importantly, we will derive the 
other equation regarding S to complete the dual equa- 
tions needed for fully describing the intensity and fre- 
quency of the microlaser. 

The coupled equations can be written in a form for 
convenient interpretation. First, the integration of Eq. 



28) using Eq. ([34]) gives 



a(t) + 7 c a(t) 



N 



2g 2 a 



4g 2 a 2 + (A + d) 2 
x sin 2 {\Jg 2 a 2 + (A + <5) 2 /4r 



(36) 



Multiplying 2a on both sides leads to the familiar rate 
equation 



n + 2j c n = — 

T 



4g 2 n 



4g 2 n + (A + 5) 2 _ 
x sin 2 {^/g 2 n + (A + 6) 2 /4r) 



(37) 



where n = a 2 (equivalent to (n) in the quantum case). 
Note the modification of the gain-loss description in 
terms of the effective detuning A' — A + 5. In the steady 
state (n = 0), n is determined by 



N 
27^ 



4g 2 n + (A + 5) 2 



x sin 2 ' g 2 n + (A + <S) 2 /4-, 



(38) 
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FIG. 12: The comparison of the quantum (thick curve) and 
the semiclassical (thin curve) calculations with gr = 0.496 
and N = 0.5 for (a) the mean photon number and (b) the 
mean frequency pulling. 



Let us turn to the integration of (^-equation, which is 

s=- Ng rm dt > 



Jo 2a(tf)r 
Ng 2 (A + 6) 
'4g 2 n + (A + 5y 



sin a/ Ag 2 n + (A + 5) 2 t 



y/4g 2 n + (A + 5) 2 i 



(39) 



Note the exact correspondence of Eq. (139)) to the quan- 
tum result, Eq. (IT9l) . The two equations, Eqs. (|38f and 
(|39| . complete the semiclassical description to correctly 
give the oscillation frequency as well as the intensity. The 
correction to the previous approach can be easily appre- 
ciated by the graphical method presented in Appendix 

m 

Since those coupled equations cannot be separated for 
individual variables, we presented the numerical calcula- 
tion of n and S in Fig. 111! They are obtained by iterating 
the integration until a converging set of (n, 8) is found. 
When compared to the quantum calculation of (n) and 
(<5„), the semiclassical calculation is in good agreement 
in the example of Fig. QT] and many others. 

Of course, the semiclassical theory is limited in the 
following aspects. First and obviously, the intrinsic sta- 
tistical properties such as p n or Ai5„ are not accessible 



by the semiclassical theory. Second, the semiclassical so- 
lution provides only the possible branch solutions if the 
multiple thresholds are involved. How those multiply seg- 
mented branches are connected cannot be traced since it 
is the quantum phenomena often modeled by the tunnel- 
ing of the photon number distribution over an effective 
potential barrier [22|. Lastly, the possible branch solu- 
tion itself may not be reliable if one deals with small 
photon numbers and strong couplings as we encounter 
in Fig. In such a regime the distinction between the 
quantum Rabi frequency g^/n + 1 and the semiclassical 
approximation g\/n becomes significant, and makes a big 
difference. 



IV. CONCLUSION 

We have formulated a quantum theory of the frequency 
pulling in the cavity-QED microlaser by considering the 
effect of the dispersive atom-photon interaction in the off- 
resonance spectra. The amount of the pulling induced by 
a small number of gain atoms turns out to be significant 
due to the strong atom-cavity coupling, which also gives 
rise to a distinctive periodic nature in the pumping de- 
pendence as contrasted to the conventional laser. The 
cavity-tuning curve exhibits a dispersion curve with no- 
table features such as the multiple-step dispersion. The 
idea of quantum frequency pulling is introduced and it 
is shown that the photon-number-dependent frequency 
pulling induces a spectral broadening, which can be a 
major source of decoherence in this laser/maser. The 
spectral lineshape deviates from a Lorentzian curve de- 
pending on the strength of the dispersive interaction in- 
volved. The present work can be used to obtain the the- 
oretical spectrum of the microlaser/micromaser over the 
entire tuning range. For the microlaser, a direct mea- 
surement of the spectrum is possible owing to the capa- 
bility of photodetection in the visible region, and such a 
measurement is expected to arrive soon in view of the re- 
cent development of the high sensitivity heterodyne spec- 
troscopy [26l428| . The concept introduced here may also 
be applicable to the other cavity-QED lasers [29l - l3l| . 

This work was supported by Korea Research Founda- 
tion (Grant Nos. WCU-R32-10045 and 20110015720). 



Appendix A: Derivation of D n and S„ 

We use 

Q n+1 = V4g 2 (n + 2) + A' 2 = n n (l + ^ 
2<? 2 



2\ 1/2 



n, 



- - Ti 



and 



1 1 2g 2 



(Al) 



(A2) 
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by assuming n much larger than unity. Then the real 
part of Eq. ([15)1 can be written as 



D„ 



,'0 n r\ . f g 2 T 

sin| — sm or 



i - 



04 



<sn,|— 1 (cosl — 1 smf — 



5 2 r 



mhi — r os vo„ 



2 7c v'n(n + 1) 



+7 c (2n + 1) 



1 — cos 



7c 
4n 



2r a sin 



2 f 9 2t \ , 7c 



\ 20„ / 4n 
while the imaginary part 5 n reads 



(A3) 



S n = -r a 
A' 



A' /0„r 

■ cos — — I sin 



n +i 



0„+i 



A' 2<7 2 A' 



O SI 3 

O n T 



x< or cos v 2 



0»T 
O n T 



, cos , 
0„ I V 2 



Q n T\ g 2 r . / O n r 



■ sin 
On \ 2 



x sin 



0„.r 



g 2 AV 

' OS 



sin(0„r) 
0„.r 



(A4) 



taking the terms up to the order of (p/O n ) 2 



Appendix B: Graphical determination of the 
steady-state intensity and frequency 

To appreciate the effect of the newly found variable 8, 
let us introduce a graphical method for the determina- 
tion of the steady-state solution (n, 8). Equations p8|) 
and (|39|) can be rewritten as the mean atom number N 
represented in (n, (5)-space as 
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FIG. 13: (color online) (a) The surface of N(n, 8) represented 
by Eq. (TETJ. (b) The surface of N(n, 8) given by Eq. flB2j. 
The positive detuning of At — 0.653 as well as gr — 0.124 
(Bl) an< 4 JcT — 0.049 are used, (c) The solution curve of (n, 8) as 
a function of N is given by the intersection of the two surfaces 
in (a) and (b). 
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and 



N = 



6 Ag 2 n+(A + S) 2 
x — ^ — x 



5 + A g' 

sin Cy/Ig^n + (A + <5) 2 r) 
^Ag 2 n + (A + 5) 2 t 



, (B2) 



respectively. From these two surfaces of ./V we can con- 
struct the evolution of (n, S) as a function of the pumping 
as in Fig. [TU] 

An example is shown in Fig. [T3] for a blue detuning 
(A > 0). The surface in Fig. [TBTa) basically represents a 
multi-branch solution of n, here depicted up to the sec- 
ond branch as a function of N. The mean photon num- 
ber or the intensity n does not change much for the range 
of 5 considered in Fig. [TBTa). in contrast to its marked 
variation with N. Thus the traditionally accepted solu- 
tion obtained by taking the line of <5 = is good enough 



for qualitative understanding of n. Determination of the 
actual nonzero S and its corresponding n requires the 
consideration of Eq. (|B2|) as in the surface in Fig.fTBTb). 
Only at the very bottom of the surface (N = 0), the 
solution reduces to 6 = line. This surface shows the 
qualitative feature of the pulling. For example, it is bent 
toward negative S as N is grown. The direction of bend- 
ing is the opposite for red detuning (A < 0), although 
not shown. 

The solution is then given by the intersection of the 
two surfaces as shown in Fig. [TBTc) . The dependence of 
n (or 5) on the pumping via the atomic flux is obtained 
by projecting the intersection onto n (or 5) plane as we 
increase N vertically. The intensity n does increase as 
expected, but the bending towards the negative S leads to 
slightly increased n compared to that of 5 = case. This 
is understandable since S reduced the effective detuning 
A' = A + 5. Owing to the periodicity in the surface given 
by Eq. ()F31|) . the solution is also periodic. 
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